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ABSTRACT 

A concise stratospheric model has been used in a Monte-Carlo analysis of 
the propagation of reaction rate uncertainties through the calculation of an ozone 
perturbation due to the addition of chlorine. Two thousand Monte-Carlo cases 
were run with 55 reaction rates being varied. Excellent conveigence was ob- 
tained in the output distributions because the model is sensitive to the uncer- 
tainties in only about 10 reactions. For a Ippbv chlorine perturbation added to 
a 1.5ppbv chlorine background, the resultant lo uncertainty on the ozone per- 
turbation is a factor of 1.69 on the high side and 1. 80 on the low side. The 
corresponding 2a factors are 2.86 and 3.23. Results are also given for the un- 
certainties, due to reaction rates, in the ambient concentrations of stratospheric 
species. 
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A. INTRODUCTION | 

M 

Ruiklel et al. (1977), hereafter referred to as Paper I, have described die | 

derivation of a concise stratospheric model designed to be used in a Monte-Carlo I 

analysis of uncertainty prop^ation. One of the greatest difficulties in carrying I 

out such an analysis is the proper characterization of the imcertainty distribution | 

of the input parameters. A one-dimensional model of the stratosphere has many i 

input parameters; reaction rates, cross sections, solar intensity, transport | 

^rameters, temperature distribution, Nj and Oi model atmosphere. The most | 

important sources of uncertainty in ozone depletion calculations, as identified | 

by the National Research Council Panel on Atmospheric Chemistry (1976) are | 

the reaction rates and transport parameters. Transport is usually represented j 

in a one-dimensional model as a diffusion coefficient. This division coefficient | 

is meant to approximate die average effect of the mean global circulation system | 

and departures from this memi. This approximation is not well enou^ under- | 

stood to allow one to quantify its uncertainty. However, the reaction rates in- | 

eluded in current stratospheric chemistry models have been the subject of many 1 

laboratory studies over the past few years and, with a few notable exceptions, | 

are rapidly becoming well enough known to make a reasonable estimate of their 1 

uncertainty possible. This paper demonstrates the application of the concise I 

model described in Paper I to the propagation of reaction rate uncertainties i 

through typical stratospheric model calculations . f 

I 

B. CALCULATIONAL TECHNIQUE I 

I 

Uncertainties in rate constants used in stratospheric models range from 
rather small (approximately a factor of 1. 1) to very large (up to a factor of 10). I 

Since all reaction rates must be positive, the appropriate probability distribu- 
tion to assume is the log normal distribution (Aitchison and Brown, 1969). The 
probability that the true reaction rate lies between k and k+Ak is given by 

p(k) = (\/2w"ok)^' exp (log k - log k^)* dk 

where log k^, the logarithm of the actual measured value, is the mean value of 

1<^ k; and a, the logarithm of the estimated uncertainty factor, is the standard 

deviation of log k. The use of this distribution is equivalent to assuming that 

log k is normally distributed. This is a physically reasonable distribution in 

that laboratory experiments typically measure log k^^ rather than k^ itself i 

(R. Watson, personal communication). 
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Figure 1 is a logical flow di^pram showing how these pr<^ahility distribu- 
tions are used in the concise stratosphere model to determine the imcertainty 
propagation properties of that model. First, the model is run to iteady-state 
(< 1 part in 5000 change in ozone concentration at any altitude peT iteration step) 
with the chosen "best set" of reaction rates. This steady-state ambient result 
includes l.Sppbv CIX resulting from Ippbv methyl chloride, (CHjCl) and 0.1 
ppbv carbon tetrachloride, (CCI 4 ) at the lower boundary added to a 0. Ippbv 
CIX lower boundary condition. A chlorine perturbation of the desired magnitude 
is then added, the model is run to steady-state, and the change in column ozone, 
alone with other parameters, is determined. A random number generator is 
then used to generate a set of random numbers e. , with zero mean and umty 
standard deviation. New values for rates kj are then determined from 

Cj = (logkj - logkoil/Oj 

When all 55 rates have been thus perturbed the model is again run to steady- 
state with 1.5 ppbv CIX. The desired chlorine perturbation is then added to tbis 
model, it is run to the new steady-state, and the resulting ozone column deple- 
tion is recorded. The values of all the species concentrations before addition of 
the perturbation are also recorded. This process is then repeated by choosing 
another set of 55 random numbers, applying each to a reaction, and re-running 
the model with and without the perturbation. When 200 cases are accumulated 
the resulting distributions for the logarithms of ozone column depletion, ozone 
column density, species densities and ratios of species densities, are fitted to 
separate normal distributions on the high and low side of the mean values. 

(This procedure recognizes the fact that, due to non-linear feedback mecha- 
nisms, uncertainty distributions of model outputs may be skewed.) Further 
groups of 200 are then run until the accumulated one sigma imcertainties on 
both sides of the distribution have converged to within one or two percent. In 
all 2000 cases were nm although satisfactory convergence for column ozone 
depletion was obtained after 1200 cases. For most of the species concentra- 
tions the convergence was significantly more rapid although for a few of the 
species in the methane oxidation chain, where reaction rates are very uncer- 
tain, only poor convergence was achieved. Table 1 shows an example of the 
cumulative lo imcertainty bound for column ozone depletion at intervals of 200 
cases. Obviously the mean would have been approached in a diffej „nt manner 
had the cases been taken in a different order. Also shown are the la uncertain- 
ties for each of the groups of 200 cases independent of all of the other groups. 
(Note that the average o for the ten 200 case grcups is greater than the o for all 
2000 cases. This is because each 200 case group must be weighted according to 
how well it approximates normality, and this weighting factor tends to be smaller 
for groups with larger values of a.) 


C. RESULTS 


Figure 2 is a hisk^niin of 2000 oases of column ozone (^letton above 
4 dlffermit altitudes, 45km, 35km, 25km, and 15km, which is the lower Imund- 
ary ctf the mmiel. Ml of these results are for a perturbation of Ijqtbv in 
the aiEyn 4 )totic limit added to a l.Spidov bacl^oimd due to CH|C1 and CCI 4 as 
described above. The horizontal axis r^resents the t^tio r of the fracticmal 
ctuuQge in ozone column with randomly chosen rates to the value chtained witit 
the "best set" of rates, plotted on a logarithmic scale. The bin size is chosen 
such that A<lc^ r) is constant. The histc^rams for ozone depletion ^bove 35 and 
45km show a slight native skewness, which is most likely an mcample of a 
saturation effect. In general, a s^uration effect can be ^ected to occur when 
the ozone destruction e^iciency is h^h because most of the chlorine is in the 
form Cl or CIO and hence is participating directly in catalytic destruction. As 
can be seen from the figure, die slight skewness at high altitudes is essentially 
gone in the column above 15 km indicating that, widiin the rai^e of depletions 
represented, the catalytic efficiency i. iu>t near the saturation value. Figure 2 
also shows distribution widths which are larger at 15 and 45 km than at the mid- 
dle altitudes. The larger width at 45km is a result of control by HOj^ catalytic 
reaotion j which are more uncertain than those of , NO^, or Cl^. This effect 
is gradually washed out as the column is extended downward. The increase in 
distribution width in the total column above 15 km is due to the uncertainty in- 
herent in the subtraction of the low altitude self-healing effect from the h^^h al- 
titude destruction and to reactions Involving HO^. 

Figure 3 shows a plot, for the case of ozone depletion above 15km, of the 
logarithm of the number of cases in each bin as a function of (log r)^ for each 
side of the distribution. In these coordinates a normal distribution with zero 
mean should give a straight line of slope l/o. The best weighted strai^t line 
fit is shown. The resulting imcertainties in column ozone depletion are shown 
In Table 1. For a Ippbv Cl^ perturbation, the lo imcertainty due to reaction 
rates edone is a factor of 1.69 on the high side and a factor of 1.80 on the low 
side. The corresponding 2o uncertainties factors are 2 . 86 and 3.23. For the 
case of continuing release of fluorocarbons 11 and 12 at their 1975 rates until 
steady-state is reached (resulting in a 5.7 H>bv Cl^ perturbation) la uncertainty 
factors of 1.7 and 1.8 on the high and low sides, respectively, were obtained. 
These correspond to 2a factors of 2 . 8 and 3. 1 and are based upon an aggregate 
of only 800 cases. 

Figures 4-8 show a series of altitude profiles of stratospheric concentra- 
tions and ratios of concentrations. For each profile a range is shown corre- 
sponding to the la confidence limits calculated from the unperturbed atmosphere 
(which has 1.5ppbv Cl^). 


3 


Figure 4 shows the lo uncertainties in O 3 and CO concentration and in CH^ 
at^ NjO mixing ratio. These latter are calculated in die model by solvii^ a 
dilfusicm equation with a fix»l lower boundary condition; thus, any uncertainty 
at 15 km has been suppress^. The uncertainty increases witti increasii^ dis- 
tance from the boundary because at higher altitudes ^e mixii^ mtio is (teter- 
mlned in lazge part by what has happ^ed below and this leads to a btdld tqi in 
uncertainty. For total o<kl nitrogen, another diffuse quantity shown in Figure 
5, the lower boundary condition is unimportant compared with the stratoqiherio 
source and the uncertainty in mixing ratio is roughly constant ^ve 15 km. The 
uncertainties in NO above 45 km aiul in HNO 3 at 15 km are constrained by the un- 
certainty in total o^ nitrc^en, since at those altitudes these species dominate 
this group. 

Uncertainties in sol^ flux and photolysis cross sections (which are not in- 
clude in the present model) have been found to make only small contributions 
to the ovexall uncertainty in chlorine induced ozone perturbations (NRC Panel, 
1976). However, this does not insure that these uncertainties will be unin^r- 
tant for other model results. Care must be taken to correctly interpret Uie un- 
certainty limits calculate for any quantity for which reaction rate uncertainty 
is small and photolysis is important. Examples of this are included in the dis- 
cussion below. 

The uncertainties in O, O 3 , O 3 column, and the O/O 3 r^io merit detailed 
examination. Ozone itself comprises virtually all of odd mcygen in the altitude 
range considered, and except at 15 km, odd oj^gen is essentially produced solely 
by O 2 photolysis (see Figure 2 of i^per I). Therefore, in tiie present model the 
uncertainty in odd oj^gen above 15 km is entirely due to the imcertainty in de- 
struction reactions. Over much of the stratosphere the principal destruction 
reaction Is NO 2 + O "*■ NO + O 2 , for which the rate is believed to be exceptionally 
wen known. At the upper and lower boundaries, reactions of greater uncertainty 
are importmit and the resulting O 3 uncertainty is foimd to be greater at tiiese 
altitudes (see Figure 4). The O/O 3 ratio only defends upon the photolysis of 
ozone and its formation reaction 0 + 02 +M-+ 03 +M for which the uncertainty 
in the rate is small. Thus, the uncertainty of this ratio is small, and taken to- 
gether with the small uncertainty in O 3 , this yields a small uncertainty in O 
(see Figure 6 ). 

Figure 7 shows the uncertainties in relative amounts of the four Cl^ species. 
In practice, the uncertainties were calculated for the species concentrations in- 
dividually, but the total Cl^ mixing ratio at any hei^t is subject to only negli- 
gible errors. This is because of (1) fixed lower boundary conditions on Cl^ and 
its source molecules, ( 2 ) the exclusion of cross section and solar flux uncertain- 
ties in the present model, and (3) the small uncertainty in O 3 colunrn density. 

The last of these is in fact the only source of uncertainty present in the model 
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which affects total Clj^. This constraint on total ejq)lains tiie small imcer- 
tainty in HCl at high and low altitudes where virtually all Cl^ is HCl. 

The Cl/C 10 ratio (see Figure 8) is dependent upon the concentrations of O, 
NO and Oj and ig>on three reaction rates (Cl-fOj -► CIO + Oj , CIO + O Cl + 

O2 , CIO + NO Cl + NOj). The relatively larger uncertainty in die NO reac- 
tion rate together with the range of NO concentration shown in Figure 5 are the 
principal cause of the large imcertainty in this ratio except at high altitudes 
where the O channel becomes the dominant path for CIO reaction and here the 
uncertainty narrows considerably. 

D. DISCUSSION 

The Monte-Carlo technique has a promising application to the analysis of 
: mospheric concentration measurements. Specifically, when the ratio of con- 
oe '.trations of two constituents is measured, and the chemical time constants 
connecting these species are short compared to transport times, a one dimen- 
sional model can be expected to predict the ratio. Specific examples are O/O3, 
CI/CIO, and NO/NOj. Figure 8 shows calculated profiles for these ratios in- 
dicating the lo confidence limits. Some care, however, must be exercised in 
directly using these results in model-measurement comparisons. The ratios, 
while chemically conti*olled, are sensitive to meteorological variability. The 
0/0 j ratio depends on a dissociation coefficient and hence on the column of 
ozone above the measurement. Cl/ClO depends O/O3 at high altitude and on 
0 , NO, and O3 at lower altitude. NO/NOj depends on local O3 . Thus, a well- 
designed experiment should measure not only one of these ratios but the concen- 
trations of toe species controlling it. Then a Monte-Carlo analysis on a re- 
stricted model which fixes toe meteorologically varying components at measured 
values will yield calculationaJ. uncertainties from which an assessment of agree- 
ment or disagreement between measurement and theory can be made. 

The information content in such contparisons depends on the relative un- 
certainty bomids which can be placed on the theoretical calculation and on the 
measurement. When the experimental uncertainties are much larger than the 
theoretical ones the information content on how well the experiment confirms 
the theory is directly proix)rtional 10 the experimental uncertainty. When the 
two imcertainty bounds are comparable maximum confirmation information is 
obtained. WTieu the experimental uncertainties are sigttificantly smaller than 
the theoretical uncertainties the possibility exists for obtaining rate information. 
However, firm conclusions from this type of information should be avoided be- 
cause of the possibility that some reaction of signific.-uice has not been Included 
in the model. 
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E. SUMMARY 


This study demonstrates the ^fect of reaction rate uncertainty on tiie re- 
sults of a one dimensional stratoi^herlc chemistry model. Since the model 
results are particularly sensitive to only about 10 of the nmre than 50 reaction 
rates included only ^bout 2^^ or -1000 Monte-Carlo cases are needed to reason- 
ably span the s^e of {mssible results. This contrasts with die more than 2^^ 
calculations claimed necessary by Stewart and Hoffert (1975). This would be 
the case only if the model were approximately equally sensitive to all of the re- 
actions included. 

The uncertainties on column ozone depletion are larger than those deter- 
mined by the NRC panel. Their r^rt gives a 2o uncertainty on ozone depletion 
due to rtiiCtion rates of 2.4 deduced from a sensitivity study of 7 key reactions. 
The Monte-Carlo technique yields approximately a factor of 3 for the 2o uncer- 
tainty due to the variation of 55 rate constants. The somewhat larger number 
obtained is not particularly surprising in view of the added reactions considered. 

This calculation provides a means of assessing the uncertainty in model re- 
sults caused by model inputs in a conqirehensive manner which includes all the 
interactions of input {diameter uncertainties. The feasibility of using Uiis tech- 
nique will always rest upon the computer time required to run the underlying 
model and upon the number of cases required to (^tain convergence in the dis- 
tributions of output values. Within this constraint, the calculation can be ex- 
tended by improving the basic model and by including uncertainties in solar flux, 
photolysis cross sections, assumed mixing ratios of constituents, ten^>erature 
and other inputs besides reaction rates. The results will ne^ U7 be redeter- 
mined if and when there are substantial changes in the rate constants or their 
uncertainties or in die basic understanding of stratospheric processes. 
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Table 1 


lo Uncertainty Limits on Column Ozone Depletion Due to a 
Ippbv Perturbation for Different Numbers of Cases 


Number of Cases 

Low Side Uncertainty 

Side Uncertainty 

200 

2.16 

1.83 

200 

1.96 

1.72 

200 

1.74 

1.90 

200 

2.20 

1.70 

200 

1.71 

1.75 

200 

1.96 

1.70 

200 

1.83 

1.85 

200 

2.08 

1.78 

200 

2.19 

1.75 

200 

1.72 

1.67 

400 

1.92 

1.75 

600 

1.86 

1.77 

800 

1.84 

1.71 

1000 

1.80 

1.69 

1200 

1.80 

1.69 

1400 

1.79 

1.71 

1600 

1.80 

1.71 

1800 

1.81 

1.71 

2000 

1.80 

1.69 
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Flgfure 3. Monte-Carlo results for ozone depletion above 15km, plotted as log n 
versus (log r)2 to demonstrate the normality of the distribution (see text). 
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Figure 5. Calculated one sigma uncertainty limits of NO, NOj and 
HNO 3 concentrations and of NO^ mixing ratio. 
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Figure 8. Calculated one sigma uncertainty limits of the ratios Cl/ClO, 
O/Oj , and NO/NO^ and of the column content. 




